Astronomy & Astrophysics manuscript no. 5628 


© ESO 2008 


February 5, 2008 





H. J. Southworth^ G. Torres'*, A.J. Penny"*, J.V. Clausen^ and D.L. Buzasi^ 



Eclipsing binaries observed witli tlie wire satellite 

I. Discovery and photometric anaiysis of the new bright AO IV eciipsing binary 

O : <A Centauri 

o 

(N 

Niels Bohr Institute, University of Copenhagen, Juliane Maries Vej 30, DK-2100 Copenhagen 0, Denmark e-mail: j vc@astro .ku.dk 

04 ' ^ School of Physics A28, University of Sydney, 2006 NSW, Australia e-mail: brunttOphysics . usyd . edu . au 

^ Department of Physics, University of Warwick, Coventry, CV4 7AL, UK e-mail: j kt@astro . keele .ac.uk 

■* Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA e-mail: gtorres@cfa.harvard.edu, 

" ^ apenny@cfa.harvard.edu 

, ^ US Air Force Academy, Department of Physics, CO, USA e-mail: Derek.Buzasi@usafa.af.mil 

in 

, Received xxx; accepted yyy 

\o ■ 

O ' ABSTRACT 

\o . 

' Context. Determinations of stellar mass and radius with realistic uncertainties at the level of 1% provide important constraints on models of 
(-H Stellar structure and evolution. 

Aims. We present a high-precision light curve of the AO IV star if/ Centauri, from the star tracker on board the wire satellite and the Solar Mass 
Ejection Imager camera on the Coriolis spacecraft. The data show that if/ Cen is an eccentric eclipsing binary system with a relatively long 
orbital period. 

Methods. The wire light curve extends over 28.7 nights and contains 41 334 observations with 2mmag point-to-point scatter The eclipse 
depths are 0.28 and 0.16 mag, and show that the two eclipsing components of if/ Cen have very different radii. As a consequence, the secondary 
eclipse is total. We find the eccentricity to be e = 0.55 with an orbital period of 38.8 days from combining the wire light curve with data taken 
over two years from the Solar Mass Ejection Imager camera. 

Results. We have fitted the light curve with ebop and have assessed the uncertainties of the resulting parameters using Monte Carlo simulations. 
The fractional radii of the stars and the inclination of the orbit have random errors of only 0.1% and 0.01°, respectively, but the systematic 
' uncertainty in these quantities may be somewhat larger We have used photometric calibrations to estimate the effective temperatures of the 
components of if/ Cen to be 10450 ± 300 and 8 800 ± 300 K indicating masses of about 3.1 and 2.0 Mq. There is evidence in the wire light 
curve for g-mode pulsations in the primary star. 
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1. Introduction 

The study of detached eclipsing binaries is of fundamental im- 
portance to stellar astronomy as a means by which we can 
measure accurately the parameters of normal stars from basic 
observational data (Andersen 1991). The masses and radii of 
the component stars of a detached eclipsing binary (dEB) can 
be measured to accuracies better than 1% (e.g. Southworth et 
al.,'2005), and the effective temperatures and luminosities can 
be obtained from spectral analysis or the use of photometric 
calibrations. An important use of these data is in the calibra- 
tion of theoretical models of stellar evolution (Pols et al. 1997 
Andersen et al. 1 19901 Ribas et al. l2000> . Comparing the prop- 
Send offprint requests to: H. Bruntt 



erties of a dEB to the predictions of theoretical models allows 
the age and metal abundance of the binary to be estimated (e.g. 
Southworth et al. l2004al and the accuracy of the predictions to 
be checked, particularly if the two components of the dEB are 
of quite different mass or evolutionary stage (e.g. Andersen et 
al. [T99Tl . 

Some of the physical effects contained in current theoret- 
ical stellar evolutionary models are poorly understood, so are 
treated using simple parameteiisations, for example convective 
core overshooting and convective efficiency (mixing length). 
Whilst these parameterisations provide a good way of includ- 
ing such physics, our incomplete knowledge of their parame- 
ter values is compromising the predictive power of the models 
(Chabover ll995l Young et al. l2001.) and must be improved. The 
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amount of convective core overshooting has been studied using 
dEBs by Andersen et al. J1990> and Ribas et al. ( 2000 1, and the 
mixing lengths in the binary Al Phoenicis has been determined 
by Ludwig & Salaris (1999 1, but these studies have not been 
able to provide definitive answers because the age and chemi- 
cal composition are usually free parameters which can be ad- 
justed when comparing the model predictions to the properties 
of the dEBs. 

Progress in the calibration of convective overshooting and 
mixing length may be made by increasing the accuracy with 
which the basic parameters of dEBs are measured, which re- 
quires much improved observational data. It is arguably more 
difficult to significantly improve the quality and quantity of 
light curves, rather than radial velocity curves, because a 
large amount of telescope time is needed for each system. 
Forthcoming space missions such as kepler (Basri et al. 2005 1 
and coROT (Baglin et al. l2001> will help to solve this problem 
by obtaining accurate and extensive light curves of a significant 
number of dEBs. 

We have begun a programme to obtain high-quality pho- 
tometry of bright dEBs with the Wide field InfraRed Explorer 
(wire) satellite, with the intention of measuring the radii of the 
component stars to high accuracy. Several of the targets in this 
program have long orbital periods Cen) or shallow eclipses 
(ARCas; Southworth et al. 2006l which has made them dif- 
ficult to observe using ground-based telescopes, ij/ Cen and 
AR Cas have secondary components which have much smaller 
masses and radii than the primary stars, so will be able to pro- 
vide strict constraints on the predictions of theoretical models. 
An additional advantage of having components with very dif- 
ferent radii is that the eclipses are nearly total, which generally 
allows the radii to be measured with increased accuracy com- 
pared to systems with partial eclipses. 

In this paper we report our discovery of eclipses in the 
bright (y = 4.1) AO IV star i/r Cen made with the wire satel- 
lite. We present a detailed analysis of the light curve show- 
ing that the orbit is eccentric and the secondary eclipse is to- 
tal. As is turns out, this is a particularly interesting system 
since its brightest component is located in the region of the 
Hertzsprung-Russell diagram between the blue edge of the 
Cepheid instability strip and the region of g-mode oscillations 
seen in slowly pulsating B stars (see, e.g., Pamyatnykh 1999 1. 
We present an analysis of the oscillations showing evidence of 
two low-frequency modes that may be interpreted as global os- 
cillation g-modes in the primary star. 

2. Photometry from space: wire, smei and hipparcos 

Since 1999 the star tracker on the wire satellite (Hacking et al. 
11999 1 has been used to observe around 250 stars which have 
apparent magnitudes V <1 (see Bruntt et al. l2005t . The aper- 
ture is 52 mm and data is acquired at a cadence of 2 Hz. Each 
star is monitored for typically 2-A weeks and some targets have 
been observed in more than one season. The primary aim has so 
far been to map the variability of stars of all spectral types, and 
recently we have done coordinated ground- and space-based 
observations. When planning the wire observations a primary 
target is chosen and in addition the on-board computer selects 



four additional bright stars within the field of view (FOV) of 
8.5° by 8.5°. The monitoring of five stars by the star tracker 
was originally designed to perform accurate attitude control. 
The rms uncertainty on the position is 0.0085 pixels or about 
0.5 arc seconds. 

When binning data to one point every 15 seconds the rms 
uncertainty per data point is 0.3 and 5 mmag for a V = and 6 
star, respectively. The filter response of wire is not well known. 
From an analysis of the light collected from stars of various 
spectral types and measured brightness we estimate the filter to 
approximate Johnson V + R. The amplitudes of variation of 
variable stars studied from the ground and with wire are in 
rough agreement with this, e.g. [Bruntt et al. (2006)| compared 
Stromgren y and wire filter observations of the same 6 Scuti 
star. 

The orbital period of wire, which has been slowly decreas- 
ing since its launch in 1999, was 93.8 minutes at the time of 
these observations, wire switches between two targets during 
each orbit in order to best avoid the illuminated face of the 
Earth. Consequently, the duty cycle is typically around 20-30% 
or about 19-28 minutes of continuous observations per orbit. 

The AO IV type-star if/ Cen was observed nearly continu- 
ously for 28.7 days from June 11 to July 9 2004 as one of 
the secondary targets for the bright giant star 6 Cen. After 
binning the light curve to a time sampling of 15 seconds we 
have 41 334 data points with a point-to-point scatter of only 
2.1 mmag. The light curve shows only one primary and one 
secondary eclipse separated in time by about 7 days, indicating 
that the orbit is eccentric and has a long period. The depths of 
the eclipses are about 0.28 mag for the primary and 0.16 mag 
for the secondary, indicating that the components have quite 
diff'erent surface fluxes. The phased light curve from wire is 
shown in Fig.^ 

The WIRE data only contain one primary and one secondary 
eclipse, so cannot be used to find the orbital period of if/ Cen. 
We therefore obtained data from the Solar Mass Ejection 
Imager (smei) camera on the Coriolis satellite, smei is an all- 
sky camera aboard the U.S. Department of Defense Space Test 
program's Coriolis satellite which was launched on January 6 
2003 into a 840 km high Sun-synchronous polar orbit. The pri- 
mary purpose of smei is to detect and track transient structures, 
such as coronal mass ejections emanating from the Sun by ob- 
serving Thomson-scattered sunlight from heliospheric plasma 
(see Webb et al. 2006 ). Jackso n et al. (20 04) provide the history 
of the SMEI design and development, and the instrument design 
and testing are described by Eyles et al. (2003) A description 
of its performance is given in Webb et al. ( l2(Xj6l . 

SMEI consists of 3 CCD cameras each with an asymmetri- 
cal aperture of 1 .76 cm^ area and viewing a 3 x 60 degree strip 
of the sky, aligned end-to-end and slightly overlapping so as 
to provide a 3 degree wide strip extending approximately 160 
degrees along a great circle with one end near the direction of 
the Sun and the other at the anti-Sun direction. The cameras 
are fixed on the zenith-nadir pointing Coriolis satellite, point- 
ing some 30 degrees above the rear Earth horizon. Thus, during 
each 102 minute orbit of the satellite the cameras are swept to 
cover about 90% of the entire sky, with a duty cycle of about 
85% from launch until the present (March 2006). The unfil- 




Fig. 1. Phased light curve of i/r Cen for a period of 38.81252 days. Data points from wire (black dots), smei (grey + symbols), 
and HippARCos (grey circles) are shown. 



tered bandpass has a response which rises above 10% between 
450 and 950 nm, with a peak response of about 50%, thus cor- 
responding roughly to a wide R photometric band. A G-type 
star with V = 10 viewed on-axis produces approximately 800 
photo-electrons during a single 4 s exposure. Each camera takes 
an image every 4 seconds (which are thus smeared roughly by 
0.2 degrees across the picture). The camera pointing nearest 
the Sun (Camera 3) is used in a mode that gives it 0.1 x 0.1 
degree binned pixels, whilst for the other two the binned pixels 
are 0.2 x 0.2 degrees. The very fast (f/1.2xf/2. 2) optics are used 
to give a point spread function (PSF), with a full-width at half- 
maximum of about 0.08 degrees but with wide wings. Whilst 
for the main project aim, the optical design is very suitable, 
for the aim of stellar photometry there are certain drawbacks, 
namely the short exposures, the optical distortion of the FOV 
and of the PSF, the large pixel size, and the variation of the 
CCD temperatures during an orbit. The main limitation, how- 
ever, is the crowding. 

We used smei data covering two years from March 21 
2003 to March 15 2005. Using only one of the three cam- 
eras (Camera 2) for these preliminary studies, the data have 
gaps varying between 87 and 140 days when the star is visi- 
ble only in one of the other two cameras. The typical accuracy 
per data point is about 16 mmag', which is significantly larger 
than the wire accuracy of 2. 1 mmag. We have 3 773 data points 
from SMEI but a significant fraction are spurious outliers (see 
Sect. 12. 1> . The grey points in the phased hght curve in Fig.Q] 
are the smei data. 

We also used 91 data points from hipparcos <ESA 19^ 
collected from February 1990 to February 1993. Whilst none 
of these data points were taken during eclipse, they have been 
useful in helping to constrain the period of ij/ Cen. 

' The SMEI data is capable of giving better accuracy. On short time 
scales an accuracy of about 1 mmag is achieved, and the analysis 
method is being refined to achieve that accuracy on longer time scales. 



2.1. Data reduction 

The WIRE dataset consists of about 1.3 million windows ex- 
tracted from the 512x512 CCD star tracker camera. Each win- 
dow is 8x8 pixels and aperture photometry is carried out. The 
basic data reduction was done using the pipeline described by 
Bruntt et al. (2005 1. 

In the SMEI dataset tfr Cen will appear on approximately 30 
out of the 20000 images from each satellite orbit. For about 15 
of these, the star will be adequately distant from the FOV edges 
(giving a total exposure in each orbit of about one minute). For 
each image used, the CCD data have bias and dark current re- 
moved. The dark flux for each pixel varies from exposure to 
exposure as the camera temperature changes. The significant 
variation in effective telescope aperture with off-axis direction 
is allowed for Glare from scattered sunlight is subtracted using 
an empirical map (Buffington, private communication, 2006). 
Further details of the glare are given in Buflington et al . (2005)| 
From the satellite pointing information each pixel in an ap- 
proximately 1 .2 X 1.2 degree area around the star is given an 
absolute sky direction, accurate to a few arc seconds. Then for 
all the binned pixels used (about 500) a PSF is fitted by a lin- 
earised least-squares method to derive the height and sloping 
base. The PSF is empirically determined from observations of 
Sirius taken when the cameras are being used in an optional 
0.05 X 0.05 degree pixel high-resolution mode. Full details 
of the data reduction are planned to be given in Penny et al. 
(2006). Further discussion of stellar photometry with smei is 
also given in .Buffington et al. (2006)| 

The smei data contain a significant number of spurious out- 
liers and for the subsequent analysis we removed these. Initially 
we used all data points to determine the approximate period. 
After phasing the data points with this period, we discarded 
data points that deviate more than 3.5 cr or 64 mmag, but only 
for data points outside the times of eclipse. During times of 
eclipse we removed only five data points that deviated more 
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Fig. 2. The wire light curve during the primary and secondary eclipses is compared to the best-fitting ebop model (dotted line). 
The residuals of the fit are shown below each panel with an increased scale. In these plots groups of five points were binned 
together. 



Table 1. wire and smei observations of ip Cen. The magnitudes 
for each source have been normalised so the outside-eclipse 
magnitude is zero. The full table comprises 44 782 datapoints 
and is available at the CDS. 



HJD - 2400000 


Magnitude 


Source 


53167.7383428 


0.000022 


WIRE 


53167.7385222 


-0.000355 


WIRE 


53167.7387016 


0.002896 


WIRE 


53167.7388847 


-0.000195 


WIRE 


53167.7390662 


-0.000541 


WIRE 


53167.7392467 


0.000129 


WIRE 


53167.7394309 


-0.001556 


WIRE 


53167.7396159 


0.000402 


WIRE 


53167.7397979 


-0.001407 


WIRE 


53167.7399805 


-0.001965 


WIRE 



than ISOmmag from the mean phased light curve. A total of 
325 spurious data points, or about 9% of the smei data, were 
discarded. 

The HippARCos photometry was extracted directly from the 
VIZIER database. The 91 data points from hipparcos have a mean 
magnitude of = 4.039 ±0.004. The mean magnitude outside 
eclipse was subtracted from each of the wire, smei and hippar- 
cos prior to subsequent analysis. 

A small part of the wire light curve is given in Table ^ 
while the complete dataset from both wire and smei is available 
in electronic format at Centre de Donnees astronomiques de 
Strasbourg (CDS). 



3. Light curve analysis 

The high quality of the wire light curve of ^ Cen has allowed 
us to derive extremely precise photometric parameters for the 
stars. The orbital period of Cen is not determined from the 
wire data, which only contains one primary and one secondary 
eclipse for this eccentric system, so we have included the smei 
data in our analysis as they contain observations during eclipses 
separated by up to two years. 

Preliminary times of minimum light were derived from the 
SMEI and wire data by fitting Gaussian functions to the eclipses 
present in the data. An orbital period of 38.8119(9) days pro- 
vides a good fit to these times of minima and the overall light 
curve of if/ Cen. The actual shape of the eclipses are however, 
quite non-Gaussian, and therefore the orbital period and the 
reference time of the central primary eclipse were included as 
fitting parameters in the light curve analysis. As seen in Table|2l 
we obtained a period of 38.81252(29) days. The phased light 
curve is shown in Fig. [2 and the details of the eclipses seen 
with WIRE are shown in Fig.|2] 

The WIRE and smei light curves were analysed with jktebop 
(Southworth et al. 2004b 2004cJ, which uses the NDE ecUps- 
ing binary model (Nelson & Davis 1972 ) as implemented in the 
EBOP code (Etzel 1981 Popper & Etzel 1981J. This model rep- 
resents the projected shapes of the components of an eclipsing 
binary as biaxial spheroids, an approximation which is easily 
adequate for ijj Cen, whose components are almost spherical. 
A significant advantage of the model is that it involves very 
few calculations, meaning that extensive error analyses can be 
undertaken quickly using standard computing equipment. For 
this application, jktebop was modified to allow the orbital pe- 
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Table 2. Photometric parameters and uncertainties for 4> Cen 
found using jktebop and Monte Carlo simulations. The first 
half of the table contains the fitted parameters and the second 
half the dependent quantities. Note that 2 453 000 days must be 
added to the time of the primary minimum, Tq. 



Parameter 




Value 


Fractional sum of the radii 


'■a + '■b 


0.065861 


± 0.000070 


Ratio of the radii 


k 


0.49737 


±0.00054 


Orbital inclination (deg) 


i 


88.955 


±0.012 


Surface brightness ratio 


J 


0.688 


±0.011 


Primary limb darkening 


"A 


0.256 


± 0.009 


Secondary limb darkening 




0.362 


±0.041 


ecoso) 


ecos w 


0.52035 


±0.00011 


esino) 


esino) 


0.19036 


±0.00097 


Orbital period (days) 


P 


38.81252 


±0.00029 


Primary minimum (HJD) 


To 


183.049310 


±0.000066 


Fractional primary radius 


r\ 


0.043984 


±0.000045 


Fractional secondary radius 


'■fi 


0.021877 


±0.000032 


Stellar light ratio 


i-B 
i-A 


0.16348 


±0.00013 


Orbital eccentricity 


e 


0.55408 


± 0.00024 


Periastron longitude (deg) 


(jj 


20.095 


±0.098 



riod and time of primary mid-minimum to be included as free 
parameters. 

Preliminary solutions of the wire and smei datasets were 
made using jktebop, fitting for the orbital ephemeris and in- 
clination, the sum and ratio of the fractional stellar radii, the 
linear limb darkening coefficients for each star, and the quanti- 
ties e cos CO and e sin o) where e is the orbital eccentricity and <jj 
is the longitude of periastron. The gravity darkening exponents 
were fixed at 1.0 and the mass ratio at 0.65; large changes in 
these values have a negligible efifect on the solution. Third light 
was fixed at zero because its value converged to a significantly 
negative (and so unphysical) value if it was included as a fit- 
ting parameter. The sizes of the residuals in the smei data are 
20 mmag. The residuals of the preliminary wire light curve fit 
were analysed for asteroseismological signatures, and the two 
strongest modes were subtracted from the data (see Sect. |5}- 
The final wire light curve (41 334 data points) has a rms resid- 
ual scatter of 2.1 mmag. Once the observational uncertainties 
had been assessed for both the wire and smei data, we were able 
to combine them to obtain the final photometric solution. The 
slight difference in the wavelength dependence of the wire and 
SMEI data is unimportant here because the diff'erence in qual- 
ity and quantity between the two datasets means that the smei 
data have a significant effect only on the derived orbital period. 
The photometric parameters found in this solution are given in 
Table 121 and the best fit to the light curve is plotted against the 
WIRE data in Fig.|2l 

3.1. Light curve solution uncertainties 

The uncertainties in the photometric parameters were assessed 
using the Monte Carlo simulation algorithm contained in jkt- 
ebop (Southworth et al. (2004bl t2004ci. In this algorithm, the 
best-fitting model light curve is evaluated at the times of the 
actual observations. Observational uncertainties are added and 
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Fig. 3. Best-fitting parameter values for the 6000 synthetic 
light curves created for the Monte Carlo analysis. 



the resulting synthetic light curve is refitted, using initial pa- 
rameter estimates which are perturbed versions of the best- 
fitting parameters of the real light curve. This is undertaken 
many times, and the uncertainties are derived from the spread 
in the best-fitting parameters for the synthetic light curves. 

The WIRE light curve was re-binned to reduce the number of 
data points by a factor of 5, and then the Monte Carlo algorithm 
was run 6 000 times. The re-binning process makes a negligible 
difference to the best fit, but greatly reduces the amount of CPU 
time required for this analysis. The reduced;^'^ of the best fit to 
the re-binned data is 1.15, indicating that there is a systematic 
contribution as well as a random contribution to the residuals 
of the fit". Inspection of the light curve in Fig. |3] suggests that 
there are slight offsets in the zero-points of some of the blocks 
of observations. The final uncertainties we quote for the light 
curve parameters found in the jktebop analysis (Table are 
the standard deviations of the spread of parameter values from 
the Monte Carlo analysis multiplied by Vl.15 to account for 
the small systematic uncertainties. 



- The equivalent number, before the subtraction of the two oscilla- 
tion modes (cf. Sect.|5), was 1.56. This indicates that it was important 
to remove these modes before the final light curve analysis 
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Fig. 13 shows the spread of the values of some best-fitting 
parameters for the synthetic light curves generated during the 
Monte Carlo analysis. The radii of the two stars are not strongly 
correlated with any other parameter or with each other, so their 
values are very reliable, but a strong correlation exists between 
eccentricity and periastron longitude; this effect has previ- 
ously been noticed in CO Lac by Wil son & Woodward (1983)| 
However, the ranges of possible values for these and the other 
parameters are still very small and the correlation between ec- 
centricity and periastron longitude has very little effect on the 
fractional radii of the stars. 

As can be seen from Table |2l the fractional radii of the 
component stars of if/ Cen are known to impressive random er- 
rors of 0.10% (primary star) and 0.15% (secondary), and all 
of the light curve parameters are extremely well determined. 
In particular, the linear limb darkening coefficient of the pri- 
mary star is known to with ±0.009, which to our knowledge is 
the most precise measurement ever in an eclipsing binary sys- 
tem. However, quoting such small random errors can be mis- 
leading because systematic errors, caused by imperfections in 
the model used to fit the light curve, may be somewhat larger 
(see below). Using the more complex logarithmic or square- 
root limb darkening laws may cause the radius measurements 
of the stars to increase about 0.2% (Lacy et al. 20Q5\. 

Deriving uncertainties using the Monte Carlo algorithm in 
JKTEBOP makes the explicit assumption that the ebop model is a 
good representation of the stars. The uncertainties in the pho- 
tometric parameters arising from this assumption can be es- 
timated by fitting the light curve with other eclipsing binary 
models, such as wink (Wood [T^71> and the Wilson-Devinney 
code (Wilson & Devinnev f PTTl Wilson 1993 1. We have at- 
tempted this and found large differences between the best fits 
for different light curve models. These differences have been 
traced to possible convergence problems with the differential 
corrections optimisation algorithms contained in the Wilson- 
Devinney code, WINK, and also the original ebop code. We have 
avoided any convergence problems by generating model light 
curves, similar to that of (/rCen, with the Wilson-Devinney 
code and using jktebop to find a best fit (jktebop uses the 
Levenberg-Marquadt minimisation algorithm as implemented 
in MRQMiN by Press et al. 1992 1. We find that the radii for the 
ebop and Wilson-Devinney models agree to within 0.1% (for 
a light curve with no eccentricity) and 0.3% (for an eccentric- 
ity of 0.5), which is comparable to the random errors given 
in Table 121 Further investigations into the true accuracy of the 
photometric parameters found from the wire light curve will be 
deferred until we have modified the we^jk and Wilson-Devinney 
codes by replacing the differential corrections subroutines with 
more modern algorithms. We will then also be able to test more 
sophisticated limb darkening laws. 

4. Radial velocities 

Surprisingly little attention has been given to the measure- 
ment of the radial velocity of such a bright star as ij/ Cen. 
The first observations appear to be those made in 1914- 
1915 at the Lick Observatory ( jCampbell 1928^ , but they show 
a large scatter even in repeated measurements of the same 



three photographic plates. The remark that the hydrogen lines 
are broad probably explains this. Four additional measure- 
ments made in 1959-1960 from Mt. Stromlo were reported 
by jBusc ombe & Morris (1961)| Those velocities also show a 
large scatter, and do not phase up particularly well with the pe- 
riod and epoch we have derived from the photometry. More 
recent measurements by Grenier et al. (1999 1 are of higher 
quality (the formal uncertainty is ~ 1 km s" ' ) and do agree well 
with the expected trend. No mention of double lines has been 
made in the literature, although from our light-curve solutions 
it is expected that the secondary will be visible in the spec- 
trum. A single-lined orbital solution based on the five Grenier 
et al. 1999 velocities (attributed here to the primary star) with 
the ephemeris, eccentricity, and longitude of periastron held 
fixed at the photometric values gives a velocity semi-amplitude 
of K - 37.4 kms"\ and a minimum companion mass of 
Ml sin i - 0.495 ±0.0 11 (Ml -i-Ma)^^-'. However, reasonable es- 
timates of the primary mass (~3.1 M©; see below) then lead to 
unreasonably small values for the secondary (~1.3 M©), which 
indicates that the semi-amplitude K is significantly underesti- 
mated, by roughly a factor of 1.4. This is most likely due to 
line blending, particularly given the large rotational broaden- 
ing reported for the star Measured values of v sin / range from 
100 kms ' to 125 kms ' (Slettebak et al. [T975] Grenier et al. 
|T999 Rover etal. UM^ . 

Absolute brightness measurements of i/r Cen 
are available in a number of photometric systems 
(Mermilliod et al. 1997> . We have collected that informa- 
tion and used the measurements in the Johnson, Stromgren, 
and Geneva systems along with colour-temperature cal- 



Gray (1992) 




Napiwotzki et al. (1993)' 




"Balona (1994)| 


Smalley & Dworetsky (1995), 


and 


.Kunzlietal. (1997)1 to 



estimate the mean effective temperature of the binary. The 
result, 10 200 K, along with our determination of the ratio J of 
the central surface brightness of the components leads to indi- 
vidual temperatures of 10450 K and 8 800 K, with estimated 
uncertainties of 300 K. These correspond to approximate 
spectral types of B9 and A2. 

We then used model isochrones from the Yonsei-Yale se- 
ries JYi et al. 20 1 Demarqu e et al. 2004^ to infer the masses 
of the components by seeking agreement with both tempera- 
tures and the ratio of the radii for the system, under the assump- 
tion that the stars are coeval. For an assumed solar metallicity 
(Z = 0.01812 in these models) we find a good fit for an age of 
290 Myr. The inferred masses are 3.1 and 2.0 Mq, with an esti- 
mated uncertainty of about 10%. These mass estimates are the 
basis for our claim above that the available velocities of i/r Cen 
give a biased value of the semi-amplitude K. 

5. Search for intrinsic variability in i// Cen 

Our estimate of the mass of the primary component of i/r Cen 
places it between the blue edge of the Cepheid instability strip 
where 6 Scuti stars are found and the region of slowly pulsating 
B-type stars (SPB) which have masses above ~ 3 Mq. To search 
for intrinsic variability in i// Cen we used the wire dataset and 
removed the data collected during the two ecUpses. We then 
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calculated the amplitude spectrum which is shown in the top 
panel in Fig.0] There is a significant increase in the noise level 
towards low frequencies. Due to the duty cycle of 30% the sig- 
nal at low frequencies is leaked to the harmonics of the orbital 
frequency of wire which is fw = 15.361 + 0.001 c/day. 

We find two significant peaks at /i = 1 .996(2) and = 
5.127(3)c/day with corresponding amplitudes of 0.23(1) and 
0.18(l)mmag, and phases 0.65(6) and 0.82(8). This provides 
the best fit to S^^jO/ sin[2n(fi ■ t + <pi)]. The zero point in time 
for the phases is Tq from Table|2l The S/N levels of the modes 
are 6.0 and 8.2, respectively. We show the amplitude spectrum 
after subtracting these two modes in the bottom panel in Fig.|?] 
The inset in the bottom panel shows an increase in the noise 
level towards lower frequencies but none of the peaks below 5 
c/day are significant. 

To improve on this and look for additional variability we 
subtracted the initial solution from the light curve, i.e. the two 
modes f\ and /2 and the fitted binary light curve. We then re- 
moved slow trends in the light curve by subtracting a sliding 
mean box with a width of 2 days (high-pass filtering). A slight 
dependence of the orbital phase and flux level was seen at the 
beginning of each orbit and we decorrelated this within the 
data. The subtracted modes were then added to the improved 
light curve. We find a significant change in the noise in the 
light curve with time. During the first week and the last week 
of observations the rms noise level is about 40% higher than 
in the middle of the run. We assigned weights to each data 
point, i, as Wtime,; 1 /o"rms' where cTnns is the rms value within 
each group of 350 data points (about three satellite orbits). 
These weights were then smoothed using a sliding box with 
a width of 1.3 days. Finally, we weighted outlier data points 
using VKoutiier,/ 0^ [1 -H |Am,|/(fl ■ (ff'Y^ with constants a - 3.5 
and b — inspired by the discussion in point 5 in Sect. 4 of 
[Stetson (1990)1 Here Am, is the differential magnitude of the 
f th data point and cr is the rms of the complete light curve. 
The final weights are Wiunn ■ VVoutiiei which are then normalized 
to unity. Using the improved light curve and the weights de- 
scribed here we find the same frequencies and amplitudes of /i 
and /2 within the uncertainties, but with the S/N improved by 
20-30%. 

The primary aim of the procedure described here was to im- 
prove the light curve by subtracting the intrinsic variation in the 
B-type component and also to improve the decorrelation with 
orbital phase before making the final light curve analysis. This 
indeed improved the;^'^ fit significantly as discussed in Sect.|3] 
The secondary aim of this procedure was to detect additional 
low amplitude modes but none were found. 

5.1. Interpretation of f\ and f2 

The two detected peaks in the amplitude spectrum coiTespond 
to periods of about 0.5 and 0.2 days. These periods could be 
present in the data due to modulation caused by the rotation 
period of the stars. To examine this possibility we must esti- 
mate the absolute radii of the component stars. To convert our 
accurate relative radii into absolute radii we estimate the semi- 
major axis from Kepler's third law with the orbital period and 



f < i 

\ \J' f \^ < 


^ ^ ii 

J W I ■^J W: 


0,1 

^ 0,0 


2 4 





5 10 15 20 25 30 

Frequency [c/day] 



Fig. 4. The amplitude spectrum of the wire light curve of ip Cen 
is shown in the top panel. Two significant peaks at fi and /2 are 
found. The aliases and the wire orbital frequency are marked. 
The bottom panel is the amplitude spectrum after subtracting 
these two modes. The inset shows the details at low frequen- 
cies; note the change in scale on the ordinate. 

the mass estimates given in Sect. |4] The result is 3.7 + 0.2 Rq 
for the primary and 1 .8 ± 0. 1 7?© for the secondary. 

We use vsin/ = 125 ± 15kms ' for both the primary and 
secondary component to find rotational periods of 1.49 + 0.26 
and 0.74 + 0.13 days, respectively. To reach a period as low as 
0.5 days the secondary star (with R/Rq - 1.8) would need to 
rotate with v sin / ^ 180 km s"' . At present we cannot entkely 
rule this out, since the spectroscopically determined rotational 
velocity is dominated by the light from the primary. 

To conclude, it seems difficult to explain the observed peaks 
in the amplitude spectrum as due to the rotation of the primary 
and possibly also the secondary. 

An alternate explanation is that the variation is caused 
by high order low degree g-mode oscillations in the primary 
star. The instability strip predicting the excitation of modes 
in slowly pulsating B-type stars (SPB) was calculated for so- 
lar metallicity by Pamyatnykh (I1999> . However, we have Tes 
- 10,450 + 300 K for the primary and this is on the cool side 
of the theoretical SPB instability strip. 

C. Aerts (private communication) kindly computed an evo- 
lution model for M/Mq = 3.1 using the new solar com- 
position ( |Asplund et al. 2004) with the cles code. The pulsa- 
tion code MAD was used to test the mode excitation (see, e.g. 
Pupret et al.2003l for more information on these codes). None 
of the modes were found to be excited for this model for Teff 
in the interval 10450 ± 300 K. Ignoring this fact we compared 
the location of the two modes /i and /a to see if we could find 
a matching evolutionary state for modes with angular degree 
Z = 1 or2. 

The modes in the pulsation model are closely spaced 
around the low frequency mode /i and we cannot find an unam- 
biguous fit. In the frequency range around the modes in the 
model are well separated but both / = 1 and 2 are possible. For 
/ - 1 there is a match of observed and computed frequencies at 
Teff = 10 250 K and for / = 2 both T^s = 10 200 K and 10 600 K 
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are possible. However we must remember that these estimates 
are for solar metallicity while the metallicity of if/ Cen is un- 
known. Accurate temperatures can be found once photometry 
of the eclipses in multiple filters has been collected and so a 
definite mode identification should be possible. Also, when the 
radial velocity curve of ^ Cen is acquired the masses should be 
constrained to better than 1 % and a detailed modelling of this 
system can be done. This should enable us to probe the extent 
of core overshoot in the component stars. 

It is important to note that the model code does not take ro- 
tation into account. |Aerts & Kole nberg ( 2005) i used the same 
evolution code to interpret observations of the photometric 
variations seen in the main sequence SPB star HD 121190 
which is about 1 500 K hotter than ij/ Cen. They stress that 
significant frequency shifts are expected for even moderately 
rotating stars due to the Coriolis force perturbing the global 
oscillation modes when the mode frequency is higher than 
half the rotation frequency. This is indeed the case for both 
modes seen in ij/ Cen. From the estimate made above we find 
f^rot = 1/^'rot - 0.67 ± 0.12c/day for the primary star when 
assuming iVot - 'orbit - 89.0°. 

Since current models do not predict oscillations at the 
Teff of the primary star in Cen another interesting pos- 
sibility for the excitation could be tidally induced modes. 
[Will ems & Aerts (2002) investigated this and found that in rel- 
atively close binary systems with highly eccentric orbits it is 
possible for / = 2 and \m\ - 2 modes to occur. For this mech- 
anism to work the pulsation frequency in the co-rotating frame 
must be an integer multiple of the binary orbital frequency. 
if/ Cen has a long period and so the orbital frequency is low, 
/orbit - 0.026 c/day. Hence, tidal interaction is most likely not 
the excitation mechanism in this system. Another explanation 
for the observed modes is excitation of retrograde mixed modes 
as suggested by [Towns end (2005) To explore these possibili- 
ties further would require detailed modelling of the ^ Cen sys- 
tem, which is beyond the scope of the current study. 

We should mention that wire data are available for the 
somewhat hotter star ( Ophiuchus (HD 152614). The spectral 
type is B8V and we find T^-s - 12 000 K from Stromgren pho- 
tometry using the calibration by Napiwotzki et al. (1993) This 
star is within the SPB instability strip predicted by Pamyatnykh 
J199 9l. The amplitude spectrum shows at least four signifi- 
cant modes at 0.82, 2.0, 4.6 and 4.9 c/day. The amplitudes are 
0.5 mmag for the low frequency mode and ^ 0.2 mmag for the 
remaining three modes. Thus, the frequency range and ampli- 
tudes are similar to ip Cen. We obtained a cles model for this 
star but the mad oscillation code does not predict any modes 
above 1.5 c/day to be excited for angular degrees I - 1 and 2. 

Lastly, we see from Fig |3 that in the range of 6 Scuti star 
variation with periods of ^1-2 hours (i.e. frequencies from 
12-24 c/day) no significant peaks are present with amplitudes 
above 0.1 mmag. In the region 14-17 c/day this claim is not 
vaUd due to the increased noise level from the aliasing. From 
the light curve solution we have Lb /La - 0.163 and so no 
6 Scuti pulsation with amplitudes above ^ 0.6 mmag is present 
in the secondary component. 



6. Conclusions 

We report the discovery that i/r Cen is a bright detached ecUps- 
ing binary. Photometry from the star tracker on the wire satel- 
lite shows that the orbit is eccentric and the period longer than 
about a month. With the photometry from the smei camera we 
have been able to determine the period of the system to be 
38.81252(29)d. 

We have fitted the wire and smei light curves using ebop 
to find precise photometric parameters for i/r Cen. Monte Carlo 
simulations show that the random error in the fractional radii of 
the stars is about 0.1%. We have attempted to assess the system- 
atic uncertainties in this modelling by obtaining solutions with 
the WINK and Wilson-Devinney codes as well as with ebop, but 
have so far been unsuccessful. Experiments have shown that 
this is probably due to convergence problems with the differen- 
tial correction algorithms in the comparison light curve codes. 
We have been able to find that the systematic uncertainty in the 
radii are at most 0.3%, and a more accurate assessment can be 
made once the comparison light curve codes have been modi- 
fied to use a different optimisation algorithm. 

Only a very limited number of radial velocity measure- 
ments for t// Cen are available in the literature, and are severely 
affected by line blending, so new high-resolution spectroscopy 
will be acquired by our group. With this data we will be able to 
measure the absolute masses and radii of the component stars 
of ip Cen with high accuracy, which should allow us to place 
strong constraints on the predictions of theoretical stellar mod- 
els for masses around 2-3 M©. This will be helped if we are 
able to determine a precise metal abundance from the spectra 
of t// Cen, but the high rotational velocity of the primary com- 
ponent may make this difficult. 

We found two significant peaks in the amplitude spectrum 
at 2.0 and 5.1 c/day with amplitudes of only 0.2 mmag. These 
frequencies are too high to be due to rotational modulation. 
We tentatively interpret them as low amplitude g-mode oscilla- 
tions in the primary star A preliminary comparison with pul- 
sation models shows that a mode identification should be pos- 
sible once Teff and the mass of primary is in hand. Thus, ifi Cen 
will be an interesting object both as a classical detached eclips- 
ing binary system and as an object for a detailed asteroseis- 
mic study. It is very interesting that our results indicate that the 
region of SPB stars as determined from theory (Pamyatnykh 
i1999j may have to be extended towards the blue edge of the 
Cepheid instability strip for main sequence stars. For further 
discussion of this point we refer to | Aerts & Kolenberg (2005)| 
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